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' Abstract 

Spin-echo experiments are often said to constitute an instant of anti-thermodynamic behavior 
in a concrete physical system that violates the second law of thermodynamics. We argue that a 
proper thermodynamic treatment of the effect should take into account the correlations between 
the spin and translational degrees of freedom of the molecules. To this end, we construct an 
entropy functional using Boltzmann macrostates that incorporates both spin and translational 
degrees of freedom. With this definition there is nothing special in the thermodynamics of 
spin echoes: dephasing corresponds to Hamiltonian evolution and leaves the entropy unchanged; 
dissipation increases the entropy. In particular, there is no phase of entropy decrease in the echo. 
, We also discuss the definition of macrostates from the underlying quantum theory and we show 

0^ ' that the decay of net magnetization provides a faithful measure of entropy change. 
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O ' 1 Introduction 



The spin echo effect [TJ has received significant attention in relation to the foundations of statistical 
mechanics [21 [3l SI O El El El Ell CD]- It arguably provides a (partial) physical realization of 
Loschmidt's velocity inversion paradox [12], albeit in the context of nuclear magnetism rather than 
of gases. For this reason, it is often stated that the spin-echo effect constitutes an instance of 
a physical system manifesting anti-thermodynamic behavior. More complex echo phenomena, for 
example, [13 1 1141 [T5] . provide a fuller realization of the Loschmidt inversion i.e., they involve the full 
many-body interactions. Moreover, echo phenomena are a testing ground for ideas on the origin 
of the irreversibility in macroscopic and mesoscopic systems. However, the issue of providing a 
quantitatively precise thermodynamic description arises even in the simplest case of Hahn echoes 

The importance of the echo effects originates from the fact that, after the inversion of spins by the 
action of external pulses, a macroscopic system seems to evolve spontaneously from a disordered 
state into an ordered one. This would appear to contradict the non-equilibrium version of the 
second law of thermodynamics. Many researchers express the opinion that this contradiction is 
only apparent; the evolution is truly irreversible and the decay of magnetization provides a measure 
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of entropy change during an echo. However, it is difficult make this statement quantitatively precise. 
The usual notions of entropy (coarse-grained Gibbs entropy or Boltzmann entropy), when applied 
to a spin system, show a decrease of entropy after spin inversion (9J [10]. Moreover, they bear no 
relation to the echo decay that is thought to provide a measure of irreversibility. In this paper, we 
emphasize the necessity of a precise thermodynamic description, where the entropy functional that 
always remains a non-decreasing function of time during an echo experiment. 

To this end, we follow Boltzmann's definition of entropy and we specify macrostates relevant 
to the system. We argue that a consistent definition of macrostates cannot separate the spin from 
the translational degrees of freedom, because they become non-trivially correlated in the course 
of a spin-echo experiment. Hence, by adapting Boltzmann's coarse-graining we define an entropy 
functional that provides a thermodynamic description of spin-echo experiments, with no violation 
of the second law. Furthermore, using simple quantum open system dynamics for the description 
of relaxation, we show that the entropy increase is indeed a monotonically decreasing function of 
magnetization decay in an echo. 

For the description of the Hahn spin echoes it suffices, as a first approximation, to treat the 
thermodynamic system as an assembly of non-interacting microsystems. Interactions are essential 
for the understanding of irreversibility, but they do not affect the issue whether there is a phase 
of decreasing entropy after spin inversion [10]. We, therefore, consider an assembly of classical 
magnetic dipoles precessing in a magnetic field. The assembly consists of N particles with magnetic 
moments rrij = m(cos#sin</>, cos#cos</>, sin#), in terms of the spherical coordinates 9 and <f>\ m 
is constant equal to gs, where g is the gyromagnetic ratio and s the magnitude of the particle's 
classical spin vector. 

The total magnetization of the system is M = £^ mj. If a constant magnetic field B is applied 
along the z axis, the equilibrium configuration at temperature /3 _1 consists of all dipoles oriented 
along the direction of the field, provided that j3mB » 1. A Impulse is then applied on the system 
rotating the magnetic moments by f so that they become oriented along the x-axis. At this moment 
(t = 0), the magnetic moment of each dipole equals rrij = m(l,0,0) and a strong magnetization 
M(0) = Nm(l,0, 0) along the x axis is measured. 

The dipoles then precess around the z axis according to the equation 



where g is the gyromagnetic ratio. Each dipole precesses with a different value Bj of the magnetic 
field, reflecting the fact that the magnetic field B is not homogeneous within the material. 
By solving Eq. ([TJ we obtain 



where u>i = gB{ is the angular frequency of precession for the dipole i. The x component of 
the magnetization is then M x {t) = m coscjjf. Assuming that the angular frequencies uii are 
randomly distributed in an interval [uj m i n ,uj max ], M x {t) rapidly becomes vanishingly small because 
of dephasing. 

At t = t, a 7r-pulse is applied on the system, so that the dipoles' configuration is transformed as 



i.e., the pulse inverts the dipoles' orientation on the x — y plane. After inversion the dipoles precess 
freely. Hence, 



rhi = guii A B, 






m(cos WjT, sin Uit, 0) — > m(cos cjj-r, — sin WjT, 0) 



(3) 




(4) 
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At t = 2t, Eq. @ predicts that rrij = m(l,0, 0), i.e., there is strong magnetization in the x 
direction, the same as at time t = 0. Apparently, the system starts from an 'ordered' state, it 
evolves into a 'disordered' one at time t = r, but after the inversion it evolves back to the initial 
'ordered' state. Hence, it seems as though the system evolves from a disordered into an ordered 
state, without any external action during the time interval [r, 2r]. 

The quantification of the latter statement requires the definition of a non- equilibrium entropy 
function for the system. If the entropy is defined in terms of the dipole degrees of freedom alone, 
then the inevitable conclusion is that the disordered state at t = r is of higher entropy than the 
ordered state at t = 2r — see Refs. [9j [10] and also Sec. 2. It follows that during the time interval 
[r, 2t] the system evolves spontaneously to states of lower entropy. This is a manifestation of 
anti-thermodynamic behavior. 

In spin-echo experiments, the cause of dephasing is the spatial inhomogeneity of the magnetic 
field. In statistical mechanics, an external magnetic field is treated as an constraint external to the 
thermodynamical system. In particular, in spin-echoes the magnetic field constraint distinguishes 
the dipoles by their position. Position is not an abstract label, it is as much a physical degree 
of freedom as spin is. The treatment the dipole degrees of freedom in isolation presupposes: (i) 
a decoupling between translation and dipole degrees of freedom, and (ii) that any correlations 
between them are insignificant. Condition (i) is not satisfied, however the translational degrees of 
freedom can be considered — to a good approximation — as a bath inducing dissipation and noise 
on the dipoles' evolution. Condition (ii) is more problematic in the sense the inhomogeneity of 
the field creates correlations between position and dipole degrees of freedom. The isolation of the 
dipole degrees of freedom is a drastic simplification because it removes all information about such 
correlations from the entropy function. 

The above argument strongly suggests that a consistent thermodynamic description of the spin- 
echo effect should involve a state space that also incorporates translational degrees of freedom. We 
demonstrate than in this case, we can define a Boltzmann entropy that accounts for the correlations 
between magnetic moments and position. Hence, we conclude that in a spin-echo experiment, the 
information of the ordered initial state is transferred into information about non-trivial correlations 
[3]. After the application of the 7r-pulse, this information in correlations is again transferred to 
information about spin order. Hence, a phase of decreasing entropy never appears. 

The structure of this paper is the following. In Sec. 2, we present the definition of an entropy 
function in terms of Boltzmann macrostates that include translational degrees of freedom; this 
definition leads to a description of spin echoes with no apparent violation of the second law of 
thermodynamics. In Sec. 3, we define this version of Boltzmann entropy function in quantum 
theory. In Sec. 4, we show that the Boltzmann entropy is an increasing function of time when 
relaxation effects are also taken into account. Finally in Sec. 5, we summarize and discuss our 
results. 

2 Boltzmann entropy and macrostates for a spin system 

In this section, we define an entropy function for spin systems that provides a consistent thermo- 
dynamic description of the spin-echo experiments. As explained in the Introduction, the key idea 
is that the entropy should also incorporate the correlations between spin and translational degrees 
of freedom, for the molecules. 

First, we examine the evolution of Boltzmann entropy for the model of Sec. 1, where the 
translational degrees of freedom are ignored. The assembly of dipoles is labeled by the abstract index 
i. Since the dipoles are independent, their statistical behavior can be described by a distribution 
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function /(#,</>), denned on the two-sphere S 2 of directions in space. 

Assuming that at t = all dipoles are oriented along the x-axis, then the distribution function 
for a single dipole approximates a delta function 632(6 — f ,</>) on the sphere S 2 . The dipole's 
precession preserves 6, hence, the delta-function evolves as 5(6 — ^ ) X]^L-oo e in ^~ Wit \ Since the 
dipoles are independent, the system's macrostate is described by the averaged density 

/(^) = l5(0-|)££e^-^). (5) 

i n 

The detailed evolution of / depends on the distribution of frequencies. We assume a Gaus- 
sian distribution around a mean frequency u> — corresponding to the mean magnetic field — with a 
deviation cr w << u>. Then, 

f(6^t) = 5(6 - |) * £ / foe-^***-^ = 5(6 - |)£ e ^)-=^. (6) 
Eq. ([6]) corresponds to a diffusive evolution equation 

df _df 2 d 2 f 



which is entropy-increasing 



For the case of dipoles that are distinguished by variables external to the system, i.e., the label 
i, we employ an averaging procedure that leads to an entropy-increasing evolution equation for the 
probability density. It follows the entropy at time t = r, when the pulse is applied, is larger than 
the entropy at time t = 0. Hence, when the system retraces its past evolution after t = r, entropy 
decreases. 

As mentioned earlier, the actual distinction of the dipoles arises from the inhomogeneity of the 
field in space, i.e., from interactions involving translational degrees of freedom of the molecules. 
Therefore, the translational degrees of freedom should also be incorporated into the definition of 
the macrostates describing the system. 

To this end, we recall Boltzmann's prescription for the macrostates and entropy of rare gases. For 
a rare gas of ./V particles, microstates correspond to points of the state space T = R 6Ar . Macrostates 
correspond to distribution functions /(x, p) on the state space of a single particle Ti = R 6 . 

The Boltzmann macrostates are constructed as follows [16]: one splits the Ti-space into cells C a 
of volume (Ax) 3 (Ap) 3 » 1. To each macrostate £ = (xi, Pi; X2, P2; • • • ; xjy, pjy) G T we assign the 
sequence of numbers n(C a ) corresponding to the number of particles such that (x, p) € C a . The 
sequence n(C a ) fully specifies a macrostate for the system. At the continuum limit, the sequence 
n(C a ) defines a probability distribution /(x, p) on T\. In what follows, it is convenient to choose 
/(x, p) normalized to unity. 

The Boltzmann entropy Sb is the logarithm of the number of microstates in each macrostate, 
hence for any sequence n(C a ), Sb ~ S a ,l nn (Co)- At the continuum limit, 

S B [f] = -N j d 3 xd 3 pf(^p) ln/(x,p) + iV- AHniV. (9) 

The key point is that Boltzmann's definition of macrostates applies to any system of ./V particles, 
even when the particles have degrees of freedom other than the translational ones. It suffices that the 
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state space of the system can be expressed as F = T^. The dipole degrees of freedom corresponds to 
a classical spin vector s of constant norm s, i.e., the corresponding state space is the sphere S 2 with 
area equal to 2s. The two-sphere is a symplectic manifold with symplectic form f2 = ^ sin 0d0 A d(j). 
Hence, the equations of motion for the dipole degrees of freedom are Hamiltonian. 

It follows that the classical state space of a single particle with spin is Ti = R 6 x S 2 . The state 
space for iV-particles with spin T = R 6iV x (S 2 ) N consists of points £ = (xi, pi, si; . . . ; xjv, pat, sjv)- 
Next, we apply Boltzmann's analysis in a straightforward way. At the continuum limit, a macrostate 
is described by a distribution /(x, p, m) on fi, normalized to unity. The Boltzmann entropy is then 
a functional of / 

S B [f] = -N J d 3 xd 3 pd 2 s /(x,p,s) ln/(x,p,s) + iV- iVlniV, (10) 

where d 2 s = sinOdOdcf). 

Next, we consider the dynamical evolution of the distribution / in accordance to the simplified 
model of Sec. 1. We assume that the particle system is subjected to an external inhomogeneous 
magnetic field B(x) = B(x)z. Here, we ignore spin-spin and spin-lattice interactions, as well as 
the action of the magnetic field on the translational degrees of freedom. We also assume that 
the translational degrees of freedom are initially in a state of thermal equilibrium. Then, the 
distribution /(x, p, s) changes in time only because of the spin's precession around the magnetic 
field's axis. In terms of the spherical coordinates 9, 4>, the time evolution law for this class of states, 
is /(x, p, 6, 4>) -> /(x, p, 6>, cj) - gB(x)), or, equivalently 

at = ~ 9B{x) W 



The key difference of Eq. (Ilip from Eq. ([6]) is that in Eq. (Ilip . x is a variable of / and not an 
external label. Hence, the necessity to average over different values of / in order to obtain a closed 
evolution equation does not arise. Substituting Eq. (fTTj) into Eq. (fTUj) we obtain 



S B [f] = T^J d'xd 3 p smOdO B(x) (J d ^^A = 0. (12) 

Therefore, we showed that the Boltzmann entropy remains constant, irrespective of whether there 
is dephasing or rephasing of the transverse magnetic moments. Moreover, the action of the -zr-pulse 
does not change the entropy. Hence, the Boltzmann entropy records no violation of the second law of 
thermodynamics during the rephasing stage of a spin-echo experiment. We therefore conclude that 
if the particle positions are included into the definition of the system's macrostates, the system's 
evolution is reversible. 

While Eq. (jlip is reversible, it does not coincide with Liouville's equation for a single dipole 
coupled to the inhomogeneous magnetic field B(x). The Liouville equation for the Hamiltonian 
H = — B(x) • m = -mB(x) sin 9 involves an additional term -mVB(x) • -g^ sin 6, in the right-hand- 
side of Eq. CLE]). ' P 

We proceed to show that the additional term is small, so that Eq. (Ilip well approximates 
the Liouville equation. Since we assume that the translational degrees of freedom are in thermal 
equilibrium, the momentum dependence of the distribution /(x, p, 9, (ft) is of the Maxwell-Boltzmann 
p 2 

type ~ e~^2M, where M is the particles' mass. Denoting by A(f) the right-hand-side term of Eq. 

(TTTj) . and by C(f) the additional term above, we find that \C(f)/A(f)\ is of the order of P^J^I M ; 
where L is the characteristic scale of inhomogeneities in the magnetic field and E the mean kinetic 
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energy of a particle. The latter is of the order of /3 _1 , hence, \C(f)/A(f)\ ~ w/3/M/L. Therefore, 
the term C(f) is negligible, if the scale of the inhomogeneities of the magnetic field is much larger 
than the thermal de Broglie wavelength XdeB ~ \J P/M of the particles. 



3 Quantum Boltzmann entropy for spin-echo systems 

The definition of the function /(x, p, s) that describes the system's macrostates incorporates Boltz- 
mann's coarse-graining on the classical state space. However, spin is fundamentally a quantum 
variable, and unlike the translational degrees of freedom it has discrete spectrum. The question 
then arises, how /(x, p, s) can be constructed in terms of the underlying quantum theory. The aim 
of this section is to provide one such construction for the distribution function. We also discuss 
some subtle points regarding the physical meaning of the corresponding macrostates. 

3.1 The description of spin-echo macrostates in quantum theory 

In derivations of the Boltzmann equation for quantum gases, the corresponding distribution function 
is usually defined in terms of the single-particle reduced density matrix |17l [18]. In particular, let 
% = £ 2 (R 3 ) (g) C 2s+1 be the Hilbert space for a single particle of spin s = ^, n = 1, 2, . . . ,. A 
system of N particles is described by vectors on the Hilbert space (®H) N . The density matrix p of 
the iV-particle system is totally symmetrized for bosonic particles and totally antisymmetrized for 
fermionic particles. The single-particle reduced density matrix on Ti is defined as p\ = Tr/^ H \N-ip. 

Given the single-particle density matrix p\, we construct different versions of the functions 
/(x, p, s) according to the theory of quantum quasi-probability distributions. In derivations of the 
quantum Boltzmann equation for gases the Wigner function is usually employed [17], mainly because 
it simplifies the calculations. A Wigner function on the single-particle state space Ti = R 6 x S 2 is 
indeed defined from the reduced density matrix p\ [19] . However, Wigner functions are not positive- 
valued in general, therefore: (i) they do not have a physical interpretation in terms of particle number 
as in Boltzmann's definition of macrostates, and (ii) they cannot be used to define Boltzmann 
entropy according to Eq. (jlOp . A coarser quasi-probability distribution that is positive-definite 
should be used instead. In general, such distributions are constructed from Positive-Operator- 
Valued Measures (POVMs) [20]. Any family of positive operators, n(x, p,s) normalized to unity 
as 



for some invariant measure dp on Ti, defines a mathematically appropriate probability distribution 



The simplest case of such a POVM is obtained from the coherent states |x, p,s) on T~L, setting 
n(x, p, s) = |x, p, s) (x, p, s| , so that 



The density /(x, p, s) is the Husimi distribution associated to the single-particle reduced density 
matrix p\ through the coherent states |x, p,s). 




(13) 




(14) 




(15) 
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The coherent states |x, p,s) are defined as the tensor product |x, p) <g) |s), where |x, p) are the 
standard coherent states on L 2 (R 3 ) and |s) := \8, </>) are the spin-coherent states on C 2s+1 |21j 

IM>= E ( S 2 + S m )cos s+ ^ 6 - S in s -^ 6 -e~ im ^ 2 \m s ). (16) 

m a =— s ^ ' 

In Eq. (|16p . |s, m s ) are the eigenstates of the S z generator in the (2s + l)-dimensional represen- 
tation of SU(2). The resolution of the unity for the spin coherent states is 

(2* + i)/^^|MKM = i. (17) 

Hence, the invariant measure d 2 s equals (2s + 1) sxn 4 and the volume of the corresponding 
two-sphere is (2s + 1). 

Since the particles are assumed independent, then the single-particle reduced density matrix 
evolves under the Hamiltonian H = —gB(x)fj,S z , where p is the particle's magnetic moment. A spin 
coherent state \9, (f>) in an external magnetic field along the z-direction evolves into another coherent 
state with parameters following the corresponding classical equations of motion, i.e., \9,<p) — > \9, <p + 
gBt). As in the classical case, the assumption that the inhomogeneity scale of the magnetic field is 
much larger than the thermal de Broglie wavelength of the particles suffices to guarantee that the 
distribution /(x, p,s) evolves under Eq. (fTT|) , 

3.2 Interpretation of the quasi-classical description 

In Sec. 3.1 we showed a suitable definition for the function /(x, p, s) from the underlying quantum 
description. However, we must elaborate here on the adequacy of the description of a quantum 
system by a classical variable. We ignore the translational degrees of freedom and focus on the 
spin ones, so that we express the POVM (|14p simply as fl(s). The classical state space for a single 
spin is the sphere S 2 . If C is a region of S 2 then we define tic = J*<7 ^ 2s n(s). The positive number 
Tr(piHc) is interpreted as an approximate probability that the spin vector lies in the region C. Thus 
the assignment C —> tic defines an approximate correspondence between operators and state-space 
regions. A necessary (but not sufficient) condition for this correspondence to be meaningful is that 
the volume [C] = f c d 2 s of the region C is much larger than unity [22]. A correspondence between 
quantum and classical observables exists only for coarse-grained phase space regions of volume much 
larger than fi. The area of the spin two-sphere is h(2s + 1), so a proper correspondence between 
operators and regions is only possible for values s >> 1. 

Therefore, it seems that it is not possible to define a classical description for particles with low 
values of spin — in particular spin s = ^. The answer to this problem lies in the remark that pi is 
not the density matrix of a single particle but the reduced density matrix in a system of iV particles. 
We express the probabilities obtained from p\ in terms of the density matrix p of the total system, 
and we obtain that Tr(pitlc) = Tr(pPc). 

p c = n c ®i®...®i + i®iic®---®i + -- - + i®i®---®nc, (is) 

is a positive operator corresponding to the proposition that "the spin of at least one particle takes 
values in the region C". TrpPc is therefore the corresponding probability. The positive operator 
P c corresponds to the region C = (Cx5 2 x...xS 2 )U(S 2 xCx...x5 2 )U...U(5 2 x5 2 x...xC) 
within the classical state space (S 2 ) N for the N spins. The volume of C is therefore 

[C] = TrP c = (N - l)(2s + 1)[C}. (19) 
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Hence, for systems with a large number of particles, the volume of the region [C] can be signifi- 
cantly larger than unity, and the approximate correspondence between positive operators and state 
space regions is meaningful. 

3.3 Quantum definition of spin-echo Boltzmann entropy. 

The distribution /(x, p,s) is positive-valued. Therefore, we can define the Boltzmann entropy as 
in Eq. CED 

S b = -nJ d 3 xd 3 p (2s + 1) 4 ^ n ^ # /(x, p, 9, <P) In /(x, p, 9, 0) + N - N In N. (20) 

We note that the expression 

Sw[h] = ~f d 3 xd* P {2s + 1) ^ ed9d ^ f(x, p, 9, <P) In /(x, P , 9, 0) (21) 

is the Wehrl entropy, associated to the coherent states |x, p, 9, (f>) [23]. Hence, 

S B = NS w [p 1 } + N - NlnN. (22) 

The description of macrostates in terms of the single-particle reduced density matrix p\ does 
not represent accurately the Boltzmann definition of macrostates, as described in Sec. 2. The 
coarse-graining corresponding to p\ is determined by the projectors Eq. (I18p . which represent the 
statement "at least one particle is characterized by values of the observables that correspond to 
the region C € Ti". In contrast, Boltzmann macrostates are defined in terms of the number of 
particles in C at a moment of time. The equivalence between these coarse-grainings requires the 
assumption that statistical fluctuations in the number of particles in C and particle correlations are 
negligible; then, the probabilities Tr{p\Pc) are indeed proportional to the number of particles in C. 
The latter distinction is not specific to the construction presented here, but it reflects the difference 
between the original derivation of Boltzmann's equation and the alternative derivation through the 
truncation of the Bogoliubov-Born-Green-Kirkwood-Yvon hierarchy of correlation functions. It is 
in fact the reason why the two approaches employ different physical conditions for the domain of 
validity of Boltzmann's equation. 

We believe that Boltzmann's coarse-graining is conceptually more satisfying, however in this 
paper we have chosen to work with the single-particle reduced density matrix. One reason is the 
significant technical difficulty in implementing Boltzmann's coarse-graining in quantum theory. In 
particular, an implementation of Boltzmann's coarse-graining from first principles requires a proof 
of decoherence (the corresponding variables behave quasi-classically) , a property that is likely to 
require significant restriction to the initial states of the system — see Refs. |22[ I24j . More impor- 
tantly, the difference above is not significant for the Hahn echoes studied here because the 7r-pulse 
inversion does not affect the part of the Hamiltonian that corresponds to the many-body interac- 
tion. Hence, a detailed description of the generation of irreversibility through spin-spin interactions 
is not necessary. 

It is important to note that the difference between the two coarse-grainings is more pronounced 
in spin systems than in gases. In principle, it is possible that in some systems the different coarse- 
grainings lead to different predictions. The reason is that in Boltzmann's coarse-graining, there is 
no kinematical restriction on the distribution function /(s) on S 2 . On the other hand, the form 
of the distribution function /(s), corresponding to the single-particle reduced density matrix, is 
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constrained by the value of the spin s. f(s) belongs to the subspaces of C 2 {S 2 ) corresponding to s. 
For example, for spin s = ^, the function f(8, <fi) = (6, 4>\pi\0, <j>) can only be of the form 

Q 

f(6, (f>) = sin 2 - + x cos 9 + r sin cos{<p + x), (23) 

for some constants x, r, x, while there is no restriction in the form of / when defined in terms 
of Boltzmann's coarse-graining. Hence, evolution equations obtained through Boltzmann coarse- 
graining may not have a representation in terms of the single-particle reduced density matrix of the 
system. This issue will be explored elsewhere, in relation to the Loschmidt echoes. 

Given the single-particle density matrix pi, one also has the option of employing the von Neu- 
mann entropy S v n(pi) = — NTr(pi log p{). In particular, the irreversible evolution equations con- 
sidered in Sec. 4 lead to an increase of the von Neumann entropy. However, given the fact that 
the Boltzmann coarse-graining in a spin system could lead to a distribution function that does not 
correspond to the single-particle density matrix, we consider that the Boltzmann entropy Eq. (|22p 
is a better candidate for the non-equilibrium entropy of the system. 



4 Entropy increase in spin echoes 

4.1 Evolution equation for the distribution function 

The description of the spin system described in Sec. 2 ignored the relaxation effects that characterize 
the evolution of nuclear spins. The inclusion of such effects requires the derivation of effective 
equations for the evolution of the distribution function /(x, p, s). 

There are two main sources of irreversibility in the evolution of the spin system: interactions 
between spins (dipole-dipole coupling) and interactions between the spin and the translational 
degrees of freedom (spin-lattice interaction). The former processes are responsible for the relaxation 
of the transverse components of the magnetization. The latter processes are responsible for the 
approach to thermal equilibrium. We assume that the translational degrees of freedom are in a state 
of thermal equilibrium, so that they essentially act as a thermal reservoir for the spin variables. 

Irreversible evolution is often described by a master equation for the single-particle reduced 
density matrix p\. To this end, one invokes a random field approximation, i.e., the assumption that 
each dipole evolves separately in a random magnetic field b(t), generated by the other particles. 
Assuming that the autocorrelation time of the random field is negligible, one employs the Markov 
approximation in order to obtain a master equation of the Lindblad type [25]. The simplest such 
master equation is [26] 

^ = -igB[S z , P ] - ^fa ([S-, [S+,p]} + [£+[&.,/>]]) 

_IZ ([§_,fiS+] - [S + ,S.p]) - ^[S z ,[S z ,p]], (24) 

where IV and Tl are phenomenological dissipation constants, corresponding respectively to the 
transverse and longitudinal components of the random magnetic field. The parameter a is deter- 
mined by the requirement that the stationary solution to Eq. (124j) is a thermal state at temperature 

/r 1 

a = e ^B_y (25) 
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Eq. (|24p yields the phenomenological magnetic Bloch equations for the macroscopic magnetiza- 
tion Mi = N{Si) 



M z 
M± 



1 , 



(2a + 1)M Z 



1, 



TigBM± - [ T T (a + -) + T L ) M±. 



(26) 
(27) 



The Bloch equations are usually expressed in terms of the spin-lattice relaxation time T\ and the 
spin-spin relaxation time T 2 , which are identified from Eqs. (|27H27p as 



Ti = [r T (2a + I)]" 1 (28) 
T 2 = [r T (a + ^) + T L }- 1 . (29) 

For the case s = ^, we obtain the solutions to Eq. (|24p 

Pll(i) = p o(to)e- r -( 1+2Q )^) + -^-(l-e- r -( 1+2 «)^)) (30) 

1 + 2a 

Poi(t) = Poi (t ) e -^(t-to)-(r T ( a +|)+r i )(^o) (31) 

poo(t) = (32) 

If we include the translation degrees of freedom and consider an inhomogeneous magnetic field, 
Eq. (|24l) generalizes to 

^fl = -igB(x)[S z ,p 1 ) - ^a(x) + [S + [S-,p!]fj 

~t" (^-'^^+] ~ {s+>s-pi\) - "^[^' [^*>pi]]> ( 33 ) 

where in this case p\ is a density matrix on the Hilbert space % = £ 2 (R 3 ) ® C 2s+1 . 

In Eq. (|33p the parameter a depends on the position x due to its dependence on the inho- 
mogeneous magnetic field B as in Eq. (|25p . However, the a term does not change the phases 
generated during time evolution, which are the important variables in the spin-echo experiment. 
Hence, assuming that the field inhomogeneities are small, we may treat a as a constant. 



4.2 Time-evolution of entropy. 

From Eq. f|33|) we construct the distribution /(x, p,s) that describes the macrostates according 
to Eq. (|15p . The corresponding Boltzmann entropy Eq. (|22p is expressed in terms of the Wehrl 
entropy. In order to compute the latter for spin s = ^ , we exploit the fact that the Wehrl entropy is 
invariant under SU(2) transformations [23j . An SU(2) transformation can bring any density matrix 
p to its diagonal form p = diag (x, 1 — x), parameterized by the single mixing parameter x £ [0, 1]. 
Then, 



S w [p] 



1 



2x - 1 



(1 - x) 2 (ln(l - x) - ~) - x 2 {\nx - \ 



(34) 



We now consider an evolution of the system according to Eq. (|33[) together with the operation 
of the two external pulses that characterize the spin-echo experiments. Initially, the system of N 



10 



dipoles occupies volume V and it is in a state of thermal equilibrium in the presence of the magnetic 
field B(x). The probability distribution factorizes as /(x, p, s) = yg(p)u t h(9,(ft), where g(p) is the 
Maxwell distribution for the momenta normalized to unity, and where Uth corresponds to a thermal 
state for the spin degrees of freedom at temperature T = j3~ l 



1 9 

uth{0, (ft) = 1 - 2q (a + sin 2 -) 



(35) 



At t = 0, a ^-pulse acts upon the system sending the positive z-axis into the positive x-axis. 
This induces a change Uth — > Uq, where 



uq(o,4>) ■- 

The system then evolves under Eq. 
evolves to 



1 + 



1 



8111 VCOS( 



(36) 



l + 2a 

The p-dependent component is unaffected, while uq 



,x 



sin + 



1 -r T (l+2a)t , 1 + (X n 

2 l + 2a K 



-^ T (l+2a)t^ 



COS ( 



+ 



1 1 



2 1 + 2a 



. e -(r T ( a +i)+r L )t sin cos r 



w«(x)t], 



(37) 



where w p (x) = gfj,B(x). At t = r, a 7r-pulse acts by inverting the spin's direction; then the system 
evolves again under Eq. ([33]) . So for t > t, 



,x 



sin - + 



+ 



2(1 + 2a) 



1 -IY(l+2a)t , 1 + a /-I _ 

2 1 + 2a 1 

e -[r T (a+i)+r L ] isin0cos r 



-r T (l+2a)t 



COS ( 



w p (x)(t-2r)]. 



(38) 



If IV > 0, Eq. (|33p does not preserve the energy, since the spin-lattice coupling transfers energy 
from the spin to the translational degrees of freedom. First, we consider the energy-preserving case 
IV = 0, and we calculate the evolution of the Boltzmann entropy Sb Eq. (jlOp . The dependence of 
Sb on Wp(x) cancels out due to the invariance of the Wehrl entropy under SU{2) transformations. 
This property is not affected by the position-dependence of the transformations. We differentiate 
Eq. ()10p with respect to time and we obtain 



S B = , - , e- r ^ [ ^ sin 2 9 cos 0(1 + In u) > 0. (39) 
2(2a + l) J 2vr r\ j - \ j 

The last step follows from the fact that the terms multiplying the positive values of cos (ft are always 
larger than the terms multiplying negative values of cos^ ( u(9,<f>) > u(9,ir — (ft), for < (ft < ir). 
Therefore, we conclude that the Boltzmann entropy is an increasing function of time, and the 
evolution is genuinely irreversible. 

Next, we consider the general case Tt > 0. We assume that the system is in contact with 
thermal reservoir of temperature /3 _1 , coupling to the translational degrees of freedom. We also 
assume that the relaxation time of the translational degrees of freedom is much shorter than T^ 1 . 
Hence, they always remain in a thermal state. The energy AE lost by the spin degrees of freedom 
is transferred to the thermal reservoir, and the entropy of the reservoir increases by an amount 
AS r > /3AE. The total entropy change AStot satisfies the inequality 

AS tot > AS B + PAE = A(S B -PH), (40) 
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Figure 1: (Color online) The minimum value of total entropy Stot — Sb — f3H plotted as an increasing 
function of time t, for Tt/{(]^B) = 0.05 and different values of temperature T. Curve (a) corresponds to 
T = 0.3.9/iB, curve (b) to T = 0.5gfiB and curve (c) to T = g^B. 

where Sb is the Boltzmann entropy and H is the mean value of the system's Hamiltonian. From 
Eq. (|3ZD we find 

^ = -T^(l-^ rT(1+2Q)< )' ( 41 ) 

where O p = y f v u p (x) = gfiB is the spatial average of the precession frequency. 

Using Eqs. (1221 l34l) . we find that Sb — f3H is an increasing function of t for all values of 
relaxation time and temperature. The behavior of Sb — /3H is plotted in Fig. 1 for different values 
of temperature. Consequently 

dStot > d(S B -/3H) > 
dt ~ dt 1 ' 

Hence, the total entropy increases in time during the spin echo experiment: there is no anti- 
thermodynamic behavior and no violation of the second law of thermodynamics. 

The above results are valid for the Markovian master equation Eq. (|38p that governs the 
evolution of the system. The predicted decay of the echo at time t = 2r, defined as 

D (t) := = e -^T(2a + l) + r L ]r = e -2r/T 2 ^ 

K 1 M x (0) V ' 

is an exponential function. 

The distribution Eq. (|37p at time t = 2r depends on r only through exponentials. Hence, 
we perform a change of variables and we express the increase AS' tot (r) of the total entropy as a 
function of the decay D(t). AStot is a strictly decreasing function of D — see Fig. 2. Hence, the 
common conjecture that the echo decay parameter provides a faithful measure of entropy increase — 
for example, Refs. [U H5] — is verified. 

4.3 Loscmhidt echoes 

Here we considered Halm spin echoes, where the external pulses invert only the evolution by the 
free spin Hamiltonian — they do not affect the spin-spin and spin-lattice interactions. More complex 
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0.0 0.2 0.4 0.6 0.8 1.0 

Figure 2: The increase in the total entropy AStot during an echo is a monotonously decreasing 
function of the magnetization decay D(r) = ^f^jp Here, curve (a) corresponds to T\jTi = 0.75, 
curve (b) to T1/T2 = 1.5 and curve (c) to T\jTi = 3. 

pulse sequences may also achieve the inversion of spin-spin couplings, whence a more complete 
realization of Loschmisdt's idea is achieved [13]. Moreover, echoes for localized excitations have 
been obtained. In general, inversion is never perfect, and the echo signal decays with time. We 
emphasize here the existence of systems with fast spin dynamics, where the decay is not determined 
by any external interactions (e.g., spin-lattice coupling) but by the reversible spin dynamics itself 
[T5] (Loschmidt echoes). In this case, the function D(t) is best described by Gaussian rather than 
exponential decay as in Eq. (f43l) . 

In this work, emphasis is given on the construction of a non-equilibrium entropy, according to 
Boltzmann, that provides a consistent thermodynamic description of spin-echo experiments. We 
showed that this entropy is a monotonic function of the echo decay D{t). We defined the entropy 
by adapting Boltzmann's coarse-graining for the rare gases into the spin context: macrostates 
correspond to a distribution function on the state space of a single particle. There seems no difficulty 
in applying the same construction to all echo experiments. In fact, we can take the analogy with 
Boltzmann's theory of rare gases one step further. We can construct an evolution equation for the 
macrostates, in analogy to Boltzmann's equation, where the spin-spin interaction is incorporated 
in non-linear collision terms. The possibility that such an evolution equation could account for the 
Gaussian decay law in Loschmidt-echo experiments is at present being explored. 

5 Discussion 

An important conclusion is drawn from the work presented here: the interpretation of the spin- 
echo experiments depends closely on the choice of macrostates one adopts for the system. In 
particular, if the macroscopic description is at the level of the spin degrees of freedom only, then the 
conclusion that the spin echo manifests an anti-thermodynamic behavior is inevitable. If, however, 
the macrostates refer also to the translational degrees of freedom, then the loss of information due 
to dephasing is compensated by non-trivial correlations between the spin and position variables. 
If the chosen macrostates accommodate such correlations, then the corresponding entropy does 
not exhibit a decreasing phase. When dissipation effects are ignored, the entropy of the system 
remains unchanged during time evolution. The evolution law for the macrostates is reversible and 
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the dephasing is generated by volume-preserving Hamiltonian dynamics. 

The main conclusion is strengthened by the consideration of relaxation effects. The usual semi- 
phenomenological Bloch equations, describing relaxation in magnetic systems, correspond to a de- 
scription of irreversible dynamics in terms of a Lindblad master equation. We showed that the total 
entropy strictly increases as a function of time under this dynamics. In effect, there is nothing ex- 
traordinary in the thermodynamic behavior of spin echoes: Hamiltonian effects such as dephasing 
do not change the entropy and the relaxation effects increase the entropy. There is no phase of 
decreasing entropy. 

Our results have some interesting implications in relation to the foundations of statistical me- 
chanics. The properties of equilibrium thermodynamics are not always reliable guides for the de- 
scription of non-equilibrium processes. In one sense, the apparent anti-thermodynamic behavior in 
the spin-echo experiments is due to the fact that the concept of entropy in equilibrium configurations 
is naively transferred into a non-equilibrium context [8]. In an equilibrium spin system, the presence 
of a net transverse magnetization is highly 'improbable', hence it can be viewed as a witness of a 
low-entropy state. If such a state were considered to arise 'spontaneously', one would then say that 
we have a violation of the second law of thermodynamics. However, in the non-equilibrium context 
there is no one-to-one correspondence between magnetization and entropy. 

Our results also suggest that macrostates are not subjective: they do not correspond to a 
description of the system in terms of variables that are accessible to us through measurements 
of macroscopic variables. In spin-echo experiments, the spin-position correlations are not directly 
accessible, and one is, therefore, tempted to ignore them in the treatment of the system. However, 
this is an approximation. It turns out that it is a drastic one: it ignores the correlations between spin 
and position degrees of freedom and, consequently, it misrepresents the thermodynamic behavior of 
the system. 

We must emphasize that our definition of the macrostates in the spin-echo system is not arbi- 
trary. Boltzmann's coarse-graining in terms of the state space of a single particle (with spin) can be 
immediately generalized to other set-ups, including relativistic systems. It is a natural definition, in 
the sense that the macrostates carry a representation of the fundamental symmetries of the system, 
and they arguably correspond to the external operations that can be effected on its constituents. 
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